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Abstract 

We present the first preliminary results on photon multiplicity measurements at LHC at forward 
pseudorapidity (2.3 < rj < 3.9) in proton-proton (pp) collisions at ^/s = 7 TeV. The multiplicity 
distribution is found to be reasonably well explained by a double Negative Binomial Distribution 
(NBD). The average photon multiplicity increases logarithmically with -^s. It is found that none 
of the models we used (PYTHIA6D6T, PHOJET, and HERWIG) could explain the data. 



1. Introduction 

One of the aims of studying the relativistic heavy-ion coUisions is to create, study 
and understand a system of deconfined states of quarks and gluons commonly known 
as the Quark- Gluon- Plasma (QGP). A set of proposed signatures provides evidences for 
the formation of such a system in these collisions. Measurements in pp collisions define a 
baseline for drawing a definite conclusion in heavy-ion collisions. Recently, measurements 
of multiplicity and pseudorapidity distributions have been reported for charged particles 
at midrapidity in pp collisions for ^/s = 0.9, 2.36, and 7 TeV by ALICE, CMS and 
ATLAS experiments at CERN-LHC [1,2]. Photon measurements provide complementary 
information. The majority of photons emitted from the reaction are decay products of 
produced particles. The multiplicity measurements at forward pseudorapidity represent 
an additional tool for testing the limiting fragmentation (LF) hypothesis which is found 
to hold at lower energies [3]. 

In this article, we report the first measurement of the multiplicity distribution of 
inclusive photons at forward pseudorapidities of 2.3 to 3.9 using the ALICE detector at 
the CERN-LHC for pp coUisions at 7 TeV. 
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2. The experiment and data analysis 

The present analysis uses data from the Photon MultipHcity Detector (PMD). The 
PMD is a preshower detector with a three radiation length thick lead converter sand- 
wiched between two planes of highly granular gas proportional counters. The information 
from one of the detector planes, placed in front of the converter is used to veto the charged 
particles, whereas the preshower data from the other detector plane is used for photon 
identification [4]. In the present analysis data from preshower plane only is used. 

We use triggered event sample (OR between the signals from the Silicon Pixel Detector 
and the VZERO, MBqr INEL [1]) which corresponds to the minimum bias inelastic 
(INEL) events. Events with vertex Z-position within ± 10 cm from the interaction point 
are analyzed. The total number of events used in the analysis after the trigger and the 
vertex selections is around 390 K. 

3. Photon reconstruction and unfolding method 

The reconstruction of photons is achieved by finding clusters of hits and then discrimi- 
nating photons and hadrons. A simple clustering algorithm with a contiguous cell search 
is employed to find the clusters. The properties of each cluster with regard to the number 
of cells (cluster Nccii) contained in it and the total energy deposition (ADC) have been 
obtained. Discrimination between photons and hadrons is made by using a threshold 
method by putting a cut-off in the cluster Ncou and ADC content. 




Multiplicity (True) 

Fig. 1. (Color online) (Left) The detector response in terms of measured and incident photon multiplicity. 
(Right-Top) The distribution of true, measured, and unfolded multiplicity of photons. (Right-Bottom) 
The ratio unfolded/true. Results are obtained using PHOJET for pp collisions at = 7 TeV for /3 = 
100. 

The method used to correct the measured raw distribution for efficiency, acceptance 
and other detector effects, is based on unfolding with x^-minimization with regulariza- 
tion, similar to what is followed for charged particles measurement [1]. In this method all 
the detector effects are described by a response matrix, A. The matrix element A^t gives 
the conditional probability that an event with true multiplicity t is measured as an event 
with the multiplicity m (see Ref. [5] and references therein for detail). Figure l(Left) 
shows the PMD detector response in terms of reconstructed photon multiplicity vs. true 
photon multiplicity within 2.3 < 77 < 3.9 inpp collisions at ^/s = 7 TeV using 1.2 M events 
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from the PHOJET event generator. The detector simulation is done using GEANT3 pack- 
age in ALICE environment. This 2-dimensional distribution forms the response matrix to 
unfold the measured distribution in data. The unfolded spectrum is found by minimizing 

where Cm is the error in the measurement g, and U is the guessed spectrum. A regular- 
ization term I3P(U) is added to the x^-function to get rid of the spurious oscillations in 
the solution. The coefficient (3 is the weight factor. Figure 1 (Right) shows the true, mea- 
sured and unfolded multiplicity distributions of photons (Top) and the ratio of unfolded 
to true distributions (Bottom) using simulated data. The ratio between the unfolded and 
the true distribution is unity within ± 10% up to a multiplicity < ~ 40 indicating 
that the performance of the method is satisfactory. 

4. Results and Discussion 

Figure 2 (Left) shows the multiplicity distribution of photons within 2.3 < rj < 3.9 
compared to the model predictions from PYTHIA6D6T, PHOJET, and HERWIG. The 
ratios between the data and the models are shown in Figure 2 (Right). It is found that 




Fig. 2. (Color online) (Left) distributions for 2.3 < ?; < 3.9, for INEL events at = 7 TeV (Solid 
magenta circles). The lines are the predictions from PYTHIA6D6T (dashed black), PHOJET (dotted 
red), and HERWIG (solid green). (Right) The ratio of the multiplicities between the data and the model 
predications at 7 TeV (PYTHIA6D6T (black triangles), PHOJET (red circles), and HERWIG (green 
squares)). The error bars for data points represent statistical uncertainties in both the figures. 

none of the models reproduces the data. At higher multiplicities {N^ > 10) all models 
under-predict the data. However, at lower multiplicities (iV^ < 10), PYTHIA under- 
predicts the data whereas PHOJET and HERWIG seems to over-predict the data. The 
multiplicity distribution of photons can be described by a double NED as shown in 
Figure 3 (Left). Studies show that a single NBD explains the data up to a multiplicity 
of ^ 30 only, at higher multiplicities {N^ > 30) it over-predicts the data. Figure 3 
(Right) shows the ^/s dependence of the average number of photons per event within 
2.3 < 7] < 3.9 for non-single diffractive (NSD) events. The data points at lower energies 
are obtained from Ref. [6] . In Ref. [6] the UA5 experiment presented the results for NSD 
events, therefore, the data points at 7 TeV, which are measured for INEL events, are 



3 



scaled to the corresponding NSD using PYTHIA simulations. The solid line is the fit to 
the data points with a logarithmic function of the form pO + pi x In ^/s. The average 




Fig. 3. (Color online) (Left) N-y distributions for 2.3 < »7 < 3.9 for INEL pp events at ^ = 7 TeV. 
The solid line shows double NBD fit. The error bars for data points represent statistical uncertainties. 
(Right) The average number of photons measured within 2.3 < r] < 3.9, for NSD events as a function of 
^/s. The solid line is a logarithmic fit to the data. 

number of photons within 2.3 < i] < 3.9 for NSD events increases with ^/s, and the 
dependence is found to be logarithmic of the form as mentioned earlier. For comparison 
our results for INEL are also shown in Figure 3 (Right). 

5. Summary 

In summary, we have presented the first measurement of multiplicity distribution of 
photons produced in pp collisions at the LHC, at forward pseudorapidity (2.3 < rj < 3.9) 
using the PMD data at ^/s = 1 TeV for INEL events. The multiplicity distribution is 
found to be reasonably well explained by a double NBD. The average photon multiplicity 
for NSD events increases with -y/s and the dependence on ^/s is found to be logarithmic of 
the form pO + pi x In -^s. The multiplicity distribution is compared with the available 
model predictions from PYTHIA6D6T, PHOJET and HERWIG and it is found that 
none of them could explain the data. 
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